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ABSTRACT 

In physics simulation chains, the PHOTOS Monte Carlo program is often used to simulate 
QED effects in decays of intermediate particles and resonances. 

The program is based on an exact multiphoton phase space. In general, the matrix 
element is obtained from iterations of a universal kernel and approximations are involved. 
To evaluate the program precision, it is necessary to formulate and implement within the 
generator the exact matrix element, which depends on the decay channel. Then, all terms 
necessary for non-leading logarithms are taken into account. In the present letter we focus 
on the decay K — > itl^vi and tests of the PHOTOS Monte Carlo program. We conclude a 
0.2% relative precision in the implementation of the hard photon matrix element into the 
emission kernel, including the case where approximations are used. 
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1 Introduction 



Semileptonic flavour changing decays, such as K — > nl fi offer a window for measurements 
of Standard Model basic couplings: quark mixing angles pQ. Moreover properties of low 
and medium energy hadronic interactions manifest themselves in such decays. It is thus 
important to keep control of decay products distributions in a form suitable for comparisons 
with data and without loosing control of the underlying quark level matrix element. 

Comparison between experimental data and theoretical predictions relies on Monte 
Carlo simulation to take into account the detector response [2J. Given today's experimen- 
tal precision, generators used in such comparisons must be based on exact phase space 
and explicit formulation of the matrix element. Otherwise the discussion of theoretical 
uncertainties in realistic applications is rather difficult. 

QED bremsstrahlung must be taken into account in these comparisons too. Infrared 
singularities cancel out in sufficiently inclusive observables. In a first approximation, the 
QED bremsstrahlung amplitude can be factorized as the Born amplitude and an emis- 
sion (eikonal or collinear) term. This can be done for a calculation at a fixed order of 
perturbation expansion, but it holds to all orders and is known under the name of expo- 
nentiation [3] (for application in Monte Carlo simulation see e.g. [HE]). In phase space 
regions where photons are collinear to charged particles, factorization theorems define the 
dominant terms of the amplitudes. This is why, bremsstrahlung can be treated to a good 
precision independently of the decay channel. As a consequence predictions, which neglect 
QED effects, represent a valid segment of the phenomenology work. 

For the decays of a given particle, matrix elements based Monte Carlo generators are 
prepared either by theorists working on effective lagrangians (also on QCD based predic- 
tions), or by experimental physicists. Since many years the PHOTOS Monte Carlo program 
[61 [T] is used for simulation of bremsstrahlung in decays. It represents a separate segment 
of the simulation chain. 

With the increased precision of new available data, such an implementation requires 
a careful discussion of its systematic errors, which has to be repeated for each decay 
mode. The phenomenological importance of approximations needs to be analyzed, and the 
factorization of QED emission terms need to be reviewed, using as a reference solutions 
based on exact phase space and matrix element for the whole decay (thus including QED 
bremsstrahlung) . 

The PHOTOS Monte Carlo program was presented for the fist time in Ref. [6]. The 
universal kernel was introduced and was shown to work, within expected accuracy, in cases 
where comparisons with first order matrix element reference simulation programs were 
available. In references [7] and [8] the notion of iteration was introduced in the PHOTOS 
program, first, for double photon emission, and later for multiple photon emission. Because 
of a rather unique order of iteration (iteration over sources providing collinear singularities 
is performed first, then construction of consecutive photon is started), the algorithm is 
compatible at the same time with an exclusive exponentiation and resummation of collinear 
logarithms. 

In reference [9] a discussion of the exact matrix element implemented within the PHOTOS 
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kernel was performed for the Z decay. The discussion was continued in Ref. [10] for the 
decay of a scalar into pair of scalars. In this study, a detailed presentation of the phase 
space parametrization was given and it was followed in [TT] by a discussion of the matrix 
elements of the W leptonic decay and the decay of virtual photons to pairs of scalars. 
Matrix element weights became available for public use with \12\ . 

Until now, discussions of matrix elements were addressing two-body decays only. In this 
case, not only phase space at Born level is particularly simple, but also photon emission is 
easy to handle. 

In the present work, we focus on the K13 decay for which a three-body kinematic is 
present already at the Born level. By studying this decay, we test not only the matrix 
element implementation effects, but also the PHOTOS phase space generation.]. 

Our study is organized as follows. In Section 2 we present the matrix elements for Ki 3 
decays, at Born level and in the case of single photon emission. Results of calculations 
based on scalar QED [T3] and on ChPT [IJJ [T5J [TB] with truncations as in Refs. [T71 [T8] 
are compared and discussed. In this context we also investigate matrix elements and their 
factorization properties. In Section 3, we consider fully differential distributions. For that 
purpose squared matrix element are reviewed and possible options resulting from physics 
assumptions (scalar QED or (TTJ [18]) are shown. We address again factorization properties, 
this time stressing features which are useful when constructing a Monte Carlo program. 
The set-up and implementations prepared for tests are presented in Section 4. Numerical 
results are collected in Section 5. Conclusions, including an estimate of PHOTOS Monte 
Carlo simulation precision for Ki 3 decays, are given in Section 6. 

2 Amplitudes for Ke3(7) decays. 

In the first sub-section we present the amplitude for the decay of charged kaons. In the 
second one the case of the neutral kaon is considered, similarities of the two cases are 
investigated as well. 

2.1 Amplitude for K~(p) — > 7r°(q) + e~(p e ) + v e (p v ) decay 
Let us start the discussion from the amplitude of the decay 



1 Parametrization of phase space in PHOTOS is explicit and exact, if the presampler for collinear emissions 
is used only along a single charge. Otherwise, starting from the moment when the photon is supposed to be 
added to a more than two-body configuration, an approximation appears when the phase space Jacobians 
for the multitude of branches are combined. This can be improved, but for test cases with second order 
matrix elements in Z decay [9] , we have found that this approximation is necessary, unless complete double 
photon emission amplitude would be installed into the program at the same time. 




(1) 



2 



taken at the Born level. With the notation of Ref. [XTf IT8] it reads 

e 2 V us F Kn (t) (p + qY 



Born 



t-M*y 

«(p e , A e )7 At (l --f 5 )v(Pu,K 



SV^Syy > - 'Mr 



(2) 



where \ e {\ v ) denote the electron (neutrino) helicity, F K7T (t) is the form factor and t = 
(p — q) 2 . Because of the relatively small mass of the K meson, t is always -C and the 
amplitude simplifies to 



Mi 



Born 



(p + qY 



V>(Pe, KhuO- - 1fs)v{pu, K) ■ (3) 

Let us now consider the amplitudes for single photon emission. We follow Ref. [11] and 




(c) (d) 

Figure 1: Feynman diagrams of the K~ — > 7r°e - z/ e 7 decay. 



we define the photon polarization states in the rest frame of the K meson. This choice 
fixes the gauge. However, our amplitudes are expressed in a Lorentz and gauge invariant 
formulation. 

Scalar QED diagrams for K~{p) — >■ ^{q)e^{p e )v e {p u ) j(k) decay are presented in 
Figure 1. The contribution of diagram (d) is proportional to 1/M^ while for other diagrams 
it is proportional to 1/M^?. For this reason, we can neglect the contribution of (d). The 
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amplitude of the K — > n°e v e ^ decay reads: 



M c = GFVus ^ {t) u(p e ,\ e ) [Q e (p + q)» 



-Q K ^]l^-ls)v(jPv,\ v ), (4) 

where Q e , Qk denote the charges of e~ and K~ respectively. To visualize its factorization 
properties this amplitude can be expressed as a sum of three gauge invariant terms: 



where 



M c = M'i + M' h + Mj n , (5) 



G F V us F Kn (t) \u[^Pe-e ~ p-e 
Mj = —(p + qY Q e - Q K 



p e ■ k p ■ k / 

tt(Pe,A e )7i«(l-T5)v(Pi/>A I/ ), (6) 

Mjj = —(p + qYu(p e ,\ e )Q 



lnO- - 1b)v(Pv, K) , (7) 



p ■ k 

u(Pe,Kh^-l5)v(Pu,K)- (8) 

The first term Mf consists of the Born-level amplitude times an eikonal factor. The second, 
Mjj, is free of soft singularity but contributes logarithmically in the collinear limit. Finally 
the third term, Mf n , is free of both soft and collinear singularities. Hence, formula (EJ 
provides a clearly structured expression of the amplitude. 

For each term fl6|7|8|) a separation into leptonic and hadronic parts is visible. It is 
encouraging, because it coincides with the structure which was useful in [11] . In our present 
work we again see that the first two parts are process independent in their emission aspect, 
and only the last non-dominant part breaks this propertjU. This is why, we expect this 
formulation to be useful for our numerical discussion. 

To prepare a comparison with amplitudes of Ref. [17] the Born and photon emission 
amplitudes (equations fl3]) and @), can be written thanks to Dirac equation as: 



M 



G FVusFKir{t) 



Born 2 



u(p e , A e ) (2 h + m e ) (1 - j 5 )v(p u , \ u ) (9) 



2 Expression ((SJ) coincides also with the similar term (JTHJ) discussed later in the amplitude of the K° 
n^l v\ decay. 
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and 



M c = u( Pe ,\ 



Q e — — ; Qk^ 



+Q f 



2p e ■ k 



p e ■ k p ■ k / 

(2 /i + m e ){l-K)v{ Pv ,\ v ). (10) 



This new formulation coincides with formula (13) of Ref. [T7|. In that paper, the form 
factor is taken at t — (see there, formulae (2) and (3)), we follow the same choice. 
The new form of the amplitude can also be splitted into two gauge invariant parts: 



where 



M c = M c r + M c n , (11) 



A/fC G F V us F Kn (t) f p e ■ e p-e 

Mj, = U(p e , A e ) Q e T - Q K 



Z \ p e ■ K p ■ k t 

(2 ^ + me )(l- l5 )v(p u ,X u ), (12) 

Mjp = U(p e , X e )Q 



2 C, ^ C 2p e -k 

(2 /i + m e )(l-j 5 )v(p v ,\ v ). (13) 

The first gauge invariant part has the form of a born-like amplitude times an eikonal 
factor. The second one does not contribute to soft singularities, but contribute to collinear 
singularities. Unfortunately formulae f lT2| [T3|) do not manifest the process independent 
form which is useful for construction of Monte Carlo programs. 

For all amplitudes presented above as well as their parts, the terms proportional to 
the electron mass were carefully kept. The results naturally extend to the case of the 
K~ — > 7r°/i~^ decay. 

2.2 Amplitude for K°(p) — > 7r + (q) + e~(p e ) + v e {p v ) decay 

The K°(p) — > 7i + (q) + e~(p e ) + v e {p u ) decay is interesting from the point of view of Monte 
Carlo discussions, even though in this case effects of bremsstrahlung are of lesser phe- 
nomenological relevance. There are not only two charged particles of different masses in 
the final state, but there is also a spectator z/ e , important from the point of view of phase 
space generation. 

The charged and neutral K decays amplitudes are quite similar. The Born level am- 
plitude reads: 

^Born = k (P + QT 

u(p e , A e )7 M (l - 75) v (p„, A„) . (14) 
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The amplitude for single photon emission is again quite short 



M° = ^^^ U^Xe) [Qeip+qY 



p e - k 2p e ■ k 

and can be expressed as a sum of three gauge invariant parts: 

M° = Mj + M° n + M° UI , 



q ■ k 



(15) 
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(17) 



(1* 



(19) 



Here Q n denotes the 7r ± charge. As in the case of K ± , only the first part (formula (|17p ) is 
infrared singular, and the third part is free of collinear singularity. 
With the help of the Dirac equation we can transform (TIB"]) into: 
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+2 



2p e ■ k 

GFV us FKnit) 

7! 



p e ■ k q ■ k 

(2 fa + m e ) (1 - 7s)v(p*, X v ) 
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u(p e , A e )7 /1 (l - 75 )u(p„, A„), 



(20) 



which coincides with formula (14) of Ref. [T7] though their form are different. If only the 
first term of our formula (1201) would be taken, 



Ml 



GfVusFktt{^) 



s/2 

4ft 
2p e ■ k 



U(p e , X e 



0.^ + 0.'" 



p e ■ k q ■ k 
(2 /i + rrie) (1 -ls)v{p v ,\ v ). 



(21) 



Then the resulting gauge invariant formula ( 12T|) is not anymore consistent with the lead- 
ing logarithm approximation for a photon emission collinear with the 7r . On the other 
hand, for the K° decay there is no collinear enhancement of photon emission along the 7T 
direction because it is not ultrarelativistic. This is why the inconsistency has no practical 
consequences for K° decay but may be of importance for the case of decays. The case 
would be then quite similar to the one of Ref. [XT] (formulas (11) and (13) there) and is 
also rather simple to fix without return to the complete formula fl20|) . Improved in that 
respect formula fl2T|) would read: 
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u(p e , A e 
Pe ■ k 



q ■ k + p e ■ k 



Qe i Qtt , 

p e - k q ■ k 



(1 - j 5 )v(p u ,X u ). 



(22) 



We investigate it, as another option for bremsstrahlung matrix element in neutral Km 
decays. 

Our formula (!20l) can be re- written also as a sum of M®, and M^,: 



M = M°, + M / ° r , 

n G FV us FK-w{t) 



k" 



y/2 \q-k + p e -k 
-Q^) u(p e , A e )7 At (l - 7 5 )u(p I /, A 



{Qnq ■ e + QePe " e) 



(23) 



where M n i is free of singularities. 



3 Differential decay probability. 

As the spin states of the Ki 3 decay products are not measurable, we concentrate on dif- 
ferential distributions obtained from amplitudes squared and summing over lepton spin 
states. 

3.1 Born and real emissions 

Let us explore squares of amplitudes given by eqs. (IT2 | [T3|) . We use the following notations: 

S = 2p e -p u ,T = 2q-p e ,U = 2q-p u . (24) 
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The Born level expression for the charged K decay reads: 
£| M c orn|2 = Gl\V us fFUt) 32 

Spin 
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For the bremsstrahlung case, the square of the amplitude is given by 

Spin 



where 

E \ M r 
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(29) 



Here ei, e 2 are two orthogonal photon polarization vectors. As expected ^2 \Mf,\ consists 
of a Born-like expression multiplied by an eikonal factor. The second and third terms, 
\ Mfj, | 2 and 2 MpMfp*, are free of soft singularities. 

Formulae for charged K~ decay and neutral K° — » 7r + e _ z/ e are similar. Obtained from 
formula (TTi 



Born I 



(30) 



El M Born! 2 = 2x^|M^ 

Spin Spin 

differs from the Born contribution of charged K decay by a factor of 2. For the bremsstrahlung 
case, the amplitude squared can be again separated into parts. The first one consistsEI of 
the squared formula ( 1221) . which reads: 



^|M°,| 2 = 2x^|M 

Spin Spin 



c|2 



+32 



G 2 F \V us \ 2 F^(t) Pe -k 



E Q 
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q ■ k + p e ■ k 



q ■ e. 
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+2 U-p„ + 



p e • k q ■ k 

Pe ■ k 



(k ■ p u (2q ■ p e + m 2 ) 



q ■ k + p e ■ k 



k ■ p v \ k ■ p e - 2q ■ kp u ■ p ( 



+ 2 Qe (Q< 

i=l,2 ^ 



P '' e '+Qj' e ' 



p e -k q ■ k 

{p v ■ €iq ■ k - q ■ €ip u ■ k)] . 



'First term can be obtained from formula (|27|28|29|) with the help of the change of variables 

p ■ t% q- e< 



Q 



K 



p ■ k q ■ k 



(32) 
(31) 
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The second part, contributing to the squared amplitude for neutral K decay, consists of 
the squared formula (123]) and its interference with formula (122]) . We do not write explicitly 
this expression here, as it is rather lengthy, and is numerically small. 

3.2 Virtual corrections 

Virtual corrections to decays can be found eg. in Ref. [17], they are of lesser impor- 
tance than the real ones and we do not recall them here. The present work is devoted to 
the discussion of real emission corrections, which are experimental condition dependent. 
Virtual corrections have to be divided into two parts. One part is large, but adds up to zero 
with real emissions in total rate thanks to Kinoshita-Lee-Nauenberg theorem. This part 
is taken into account by the PHOTOS Monte Carlo code. The other part has to be included 
in the form-factor and incorporated to the Born level matrix element. For this purpose, 
real emission amplitude squared need to be integrated over photon momentum to control 
the sum rule to the level of complete first order. This can be performed in an approximate 
wa}H (as in Ref. [19]) or in an exact manner, following phase space parametrization as 
used in PHOTOS. Integration (analytic or numerical) need to be performed. This solution 
has to be adopted, once experimental precision approaches - precision level. Still another 
type of solution, non exploiting the sum rule of Konoshita-Lee-Nauenberg can be useful. 
We mention such possibility in footnote 



4 Monte Carlo Simulation 

We use the TAUOLA [20] code to generate Born level K decay samples. For practical reasons 
this generator is suitable for our purpose once the r decay matrix element is replaced with 
the one of the K L3 decay, and appropriate adjustment of masses and particle identifiers 
are performed. Semileptonic decays of r's are suitable for such adaptation^. Another 
advantage is that we can then guarantee full control of parameters used in our Born level 
generator and in matrix element of PHOTOS correcting weights. 

We gain, because PHOTOS is ready to use with TAUOLA. The two programs share technical 
elements. This is convenient when non-factorizable parts of matrix elements are installed. 
The interface to HEPEVT event record of the two programs offers an easy access to our testing 
tool MC-TESTER [2TJ [22] . The prepared plots are then compatible with our previous studies. 
From the user point of view, numerical results collected for this paper, and available also 
in graphic form from the web-page [22], can be of interest for benchmarking K decay 
generators independently, whether they are coded in FORTRAN or in C++. 

4 But then, no gain beyond Konoshita-Lee-Nauenberg can be achieved. In K e z decay there is no Coulomb 
effect. 

5 TAUOLA semileptonic decay channel offer an alternative crude phase space generator for radiative 
corrections, which may become useful in the future, especially if complete virtual corrections are to be 
included and studied. 
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Figure 2: Distributions of invariant masses, in GeV (GeV/c 2 , c = 1) units, constructed 
from the products of the decay K~ — > 7r°e~P e , at Born level. As in [21] V us = 0.2252 is 
used, we take however F K7T (t) = 1. This is an acceptable approximation for our purposes 
and consistent with Ref.|17j. 
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The squares of Born level amplitudes (!25|) for decays K~ — » 7r°e~i/ e , and .fT — > 7i + e~u e 
are rather easy to implement into TAUOLA. The numerical results, as histograms of invariant 
masses constructed from pairs of final state decay products, are shown in Fig. [2j We have 
compared these results with the one of Rcf. [25J and reasonable agreement was found. This 
comparison is sufficient for tests, because the matrix element is rather simple and TAUOLA 
itself is well tested. Note that detector acceptance effects were included in Ref. [25], 
therefore the corresponding figures do not coincide in all details with our Fig. [2j 

5 Numerical results 

In all tests presented in this paper we use samples of 100 milion events. We refer to 
standard PHOTOS whenever we use its publicly available FORTRAN version 2.15, or any other 
version which yields equivalent results. In particular, identical results are available (as 
default option) from the C++ PHOTOS [26J, version 3.0 or higher. One of the goals of the 
present work is to provide a systematic error for these widely used versions. 

As a first step, we perform a technical test. For the decay of charged K we have 
compared results of the standard PHOTOS with the PHOTOS version of [10J. It was the first 
version where the multiphoton phase space generation for final state of a single charged 
(and scalar) particle was exact, which was provided with the help of explicit phase space 
Jacobians. It is publicly available starting from PH0T0S++, version 3.3. We could see that 
the numerical effect is small. As expected, the difference is below 0.05 % if calculated with 
respect to the total rate. We obtained differences at the level of 10 % in regions of phase 
space contributing at the level of 10~ 4 to the total rate. 

In a second step, we have checked the contribution from collinear photon emission 
region. To this end we have selected only those events, used in the previous comparison, 
where the photon-electron pair invariant mass was at most 0.01 of their energies product 
(taken in the rest frame of K). We obtained perfect agreement, with no statistically 
significant differences. 

Only then, results presented in our article were prepared. Throughout the paper, we 
use our testing program MC-TESTER [2D 122] • The two colored (grey) lines correspond to the 
compared generation samples^]. To define the boundary of the real emission phase space 
we use the photon energy in the rest frame of the decaying kaon: it is set at 0.005 of the 
decaying kaon mass. 

We have repeated the same comparison as discussed previously but when the complete 
scalar QED matrix element is installed. As we can see from Fig. [3] the numerical effect of 
the Matrix Element and exact phase space implementation is rather small, visible only at 
the ends of the spectra (contributing at the level of 10~ 3 to the total rate) where relative 

6 The distributions of Lorentz invariants constructed from outgoing particles are shown. Additional 
information is available on plots of the web-page [23] . The SDP (shape difference parameter) obtained from 
MC-TESTER represents an exclusive surface (normalized to unity) under the green and red distributions. 
The statistical error (calculated in a rather simplified way) is subtracted from this difference. The black 
line is the ratio of the distributions. 
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differences are sizable. There, matrix element effects are rather substantial and should be 
kept in mind in some contexts, e.g. for generating background to other decays. 





Figure 3: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, con- 
structed from the decay products in K~ — > 7r°e~z/ e channel. The most sensitive invariants 
to photon energy are plotted. The red (darker grey) line is standard PHOTOS, the green 
is with exact Matrix Element. The fraction of accepted bremsstrahlung events is (7.371 
±0.0027) % in standard PHOTOS run and (7.4127 ± 0.0027) % when the matrix element is 
used. 

For the case of K~ — > 7r°/x — logarithmic corrections are of course much smaller. That 
is why non leading terms contribute to photon spectra in a relatively larger manner, and 
seemingly much larger effects can be seen in Fig. H] than in Fig. |3j Nonetheless it is not 
more significant numerically. Only 0.45 % of events enter the plots for the muonic channel; 
twenty times less than in the electron case. We can conclude that standard PHOTOS works 
sufficiently well for the K~ — > 7r°/i~z/^ decay, if one is interested in 0.2 % precision limits. 

The K° — > 7r + /~z// case is technically more interesting as two kinematical branches are 
present in the crude level phase space generator. 

In Fig. [5] we compare standard PHOTOS with a version where the scalar QED matrix 
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element is installed. The approximation in the phase space is still present. Only in Fig. [6] 
we use single channel presampler for the phase space generation and the phase space is 
exact. The effect of phase space Jacobian approximation is rather small. 

By comparing Figs. [7] and [6] we can see that the standard PHOTOS is much closer to the 
result of scalar QED matrix element than to one of (12 ip . The bulk of the difference is due 
to the non-compatibility of formula (121]) with collinear logarithms due to emission from 
charged pion (we can see it by comparing Fig. [8]). This is of course beyond the framework 
of approximation at use, but it is nonetheless of some interest, to understand the origin of 
the residual differences between possible options for the matrix element and how to fix it 
in gauge invariant way, but without study of the whole matrix element. 




Figure 4: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, 
constructed from the decay products in K~ — > n^^v^ channel. Invariants most sensitive 
to photon energy are shown. The red (darker grey) line is standard PHOTOS, the green 
is with exact Matrix Element. One could conclude that the effect of matrix element 
introduction is not small in this case. However, only a small fraction of events enter this 
plot (0.4113 ± 0.0006) % for standard PHOTOS and (0.4445 ± 0.0007) % for the case with 
matrix element. The difference is well below 0.1 % when compared to the total rate. The 
two distributions coincide in the soft photon region. 
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Figure 5: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, con- 
structed from the decay products in K° — > 7r + e~z/ e channel. The most sensitive invariants 
to photon energy are plotted. Two kinematical branches are used, thus the phase space 
is not exact. The red (darker grey) line is standard PHOTOS, the green is with exact Ma- 
trix Element. The fraction of accepted bremsstrahlung events is (8.6398 ± 0.0029) % in 
standard PHOTOS run and (8.6913 ± 0.0029) % when the matrix element is used. 
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Figure 6: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, con- 
structed from the decay products in K° — > 7r + e~z/ e channel. The most sensitive invariants 
to photon energy are plotted. Single kinematical branch is used, thus the phase space is 
exact. The red (darker grey) line is standard PHOTOS, the green is with exact Matrix Ele- 
ment. The fraction of accepted bremsstrahlung events is (8.6398 ± 0.0029) % in standard 
PHOTOS run and (8.6958 ± 0.0029) % when the matrix element is used. 
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Figure 7: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, con- 
structed from the decay products in K° — > 7r + e~z/ e channel. The most sensitive invariants 
to photon energy are plotted. Single kinematical branch is used, thus the phase space is 
exact. The red (darker grey) line is standard PHOTOS, the green is with exact Matrix Ele- 
ment. The fraction of accepted bremsstrahlung events is (8.6398 ± 0.0029) % in standard 
PHOTOS run and (8.5235 ± 0.0029) % when matrix element ([21]) is used. 
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Figure 8: Distributions of scalar Lorentz invariants, in GeV 2 (GeV 2 /c 4 , c = 1) units, con- 
structed from the decay products in K° — > 7r + e~z/ e channel. The most sensitive invariants 
to photon energy are plotted. Single kinematical branch is used, thus the phase space 
is exact. The red (darker grey) line is standard PHOTOS, the green is with exact Matrix 
Element. The fraction of accepted bremsstrahlung events is (8.6398 ± 0.0029) % in stan- 
dard PHOTOS run and (8.5928 ± 0.0029) % when matrix element ( l2Tj) with improvement of 
formula ( 122|) is used. 
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6 Summary 



In this study we have adopted the matrix element of p2] for the emission kernel of PHOTOS 
Monte Carlo [7J El |10] . We have investigated semileptonic decays of kaons into 7r and lepton 
neutrino pair. After modification for these channels, PHOTOS features exact matrix element 
(three options) and exact (or default) phase space. We have evaluated the numerical size 
of the missing terms in publicly available version of PHOTOS. The difference is of the order 
of 0.2 %, thus rather small, except for the distribution of lepton photon pair invariant mass 
spectrum, where the difference is sizable at the high end of the spectrum in case of K° 
decay. 

On the technical side, we have also checked the algorithm. We have compared the case 
when the pre-sampler is active for possible collinear singularity along lepton only (for which 
the phase space is exact and phase space Jacobians are explicit), and the case when both 
pre-samples for singularities along lepton and charged 7r directions are active. Such studies 
for more than 2 body decays were not documented until now. We found the differences to 
be below 0.05 %. Our work is supplemented with a larger set of figures, which are available 
from the web-page |23j. 

We conclude that for PHOTOS version 2.15 or higher, and for Ki 3 decays, the precision 
level with respect to matrix elements based simulation is 0.2% or better. This conclusion 
extends naturally to the case of multiple photon emission. We have identified the factor- 
ization properties of the matrix element, which were the reason why the differences were 
small. On the other hand, our error estimation is not complete. It does not include any 
discussion of uncertainty in the matrix elements due to assumptions of the models used 
for their calculation. Also, non leading effects of virtual corrections, which are expected 
to contribute at a similar 0.2 % level, are not discussed. They depend on the details of 
hadronic interactions and can not be tackled in discussion of bremsstrahlung only. 
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